Study of Chirality in the Two-Dimensional XY Spin Glass 
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We study the chirality in the Villain form of the XY spin glass in two-dimensions by Monte 
Carlo simulations. We calculate the chiral-glass correlation length exponent i/cg and find that 
UcG = 1.8 ± 0.3 in reasonable agreement with earlier studies. This indicates that the chiral and 
phase variables are decoupled on long length scales and diverge as T — > with different exponents, 
since the spin-glass correlation length exponent was found, in earlier studies, to be about 1.0. 
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Ever since VillaiJ pointed out the existence of a dis- 
crete reflection symmetry (in addition to the ordinary ro- 
tation) in frustrated vector spin systems, there has been 
considerable interest in the similarity or difference be- 
tween the behavior the variables corresponding to these 
two symmetries - chiralities (which are quenched in vor- 
tices), corresponding to reflection, and spins, correspond- 
ing to rotation. In part this interest comes from the 
observation that spin glasses seem to be in the Ising uni- 
versality class, though many of them should be described 
quite well by a Heisenberg model. This has led to spec- 
ulations that chiralities and spins order differently and 
that the Ising behavior seen in experiments might indi- 
cate the existence of a chiral glass phase in the absence of 
spin glass orderingnL|-jrhis view is supported by some nu- 
merical simulationsou. However, in spite of many studies, 
the problem still remains controversial. 

Kawamura and TanemuraEI studied the two dimen- 
sional XY spin glass by a domain wall renormalization 
group technique and were the first to present evidence 
that the chiral and spin glass correlation length expo- 
nents are different in two-dimensions. They also reported 
Monte Carlo simulations whiph supported their claimH. 
A little later, Ray and Mooreo also reported Monte Carlo 
results which indicated that in two-dimensions the chiral 
and spin glass correlation length exponents were indeed 
different. They estimated Vsg — 1-0 and i>cg — 2.0, 
the former result beina-|in good agreement with earlier 
work of Jain and Young3. More recently, KawamuraQ re- 
ported Monte Carlo simulations of the three dimensional 
XY spin glass and claimed that there is indeed a stable 
chiral glass phase. This would appear to have settled 
the issue. However, recent analytic work on the one- 
dimensional ladder latticeEl and on the two-dimensional 
system with a special choice of disorderB points in the 
opposite direction. Given this controversy it seems a rea- 
sonable time to study the two-dimensional system once 
again numerically. 

The earlier work used a representation of the model in 
terms of the phases of the XY model, see Eq. (|^) below, 
and vortices are expressed in terms of correlations of the 
phases around an elementary square. However, vortices 
are only well defined when the nearest-neighbor spin-spin 
correlation function is large, which means that the tem- 
perature is already quite low. It is therefore difficult to 
study vortex correlations over a large range of temper- 



ature. In our work, we use a different representation 
of the model, expressed directly in terms of the vortices 
themselves. As a result it is possible to study the vor- 
tex correlations over a larger range of temperature than 
before. 

The model used in the simulations is an XY spin glass 
in which the interactions, J^', have values ±J. The 
Hamiltonian is 
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where the Si are two component vectors of unit length. 
This can also be written as 
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where is the angle (phase) the XY spin makes with 
a fixed direction, and the Aw take values and tt with 
equal probability, corresponding to Jw = 1 and —1 re- 
spectively. We take the sites, I, to lie on a square lattice 
of size L X L, and the interaction is between all nearest 
neighbor pairs, counted once. 

As discussed above, it is easier to study the model in 
terms of vortices. To do so we first replace the cosine in 
Eq. (0) by the Villain periodic Gaussian function, i.e. 
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in units where Boltzmarmjls-constant is unity. Performing 
standard manipulationsallJ one finds that the partition 
function of the Hamiltonian in Eq. (^ is the same (apart 
from an unimportant smoothly varying prefactor) as that 
of the following Hamiltonian, 
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where the vortices {n^} run over all integer values, sub- 
ject to the "charge neutrality" constraint 
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and G{i — j) is the vortex interaction, 
G{i — j) J M--^ 1 — exp[ik • (ri 
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Note that the Fourier transform of the vortex interaction 
is ~ k~^ for smaU k which corresponds to a long range 
logarithmic interaction in real space. The vortices sit on 
the sites i of the dual lattice which are in the centers of 
the squares of the original lattice. The bi are given by 
(l/27r) times the directed sum of the quenched random 
variables, Aui , on the links of the original lattice which 
surround the site i of the dual lattice. They satisfy a 
constraint similar to Eq. (^, 
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replicas of the system in parallel with the same realiza- 
tion of the disorder and compute the overlap between the 
states in the two replicas, 



i=l 



(11) 



Here, and in the rest of the paper, the sum over sites i 
and j on the dual lattice is only over those sites occupied 
by a chirality. The standard spin glass order parameter 
is just [{q}], where the angular brackets denote the ther- 
mal averages and the square brackets the average over 
disorder. Two useful quantities are the Binder moment 
ratio gcG and the chiral glass susceptibility Xca defined 
by 
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and 



From now on we set the interaction strength, J, to be 
unity. 

For Aw = or TT there are two kinds of plaquettes on 
the original lattice, or equivalently two kinds of sites on 
the dual lattice: there are unfrustrated sites, on which bi 
is integer, and frustrated sites, on which bi is half-integer. 
If Aui has equal probability to takes values and tt then 
half the sites on the dual lattice will be frustrated and 
half unfrustrated on average. In the ground state, most 
unfrustrated sites will have rii = 0, while most frustrated 
sites will have rii = ±1/2. These Ising-like variables are 
precisely the chirality variables of Villain. It seems rea- 
sonable that the essential physics will be preserved if one 
just keeps these chirality variables, i.e. one fixes Ui to 
be zero on the unfrustrated sites and only allows rii to 
be ±1/2 on the frustrated sites. Thus one has an Ising 
model with a long range antiferromagnetic coupling, 
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where, for convenience, we represent the chiralities by 
Ising spins, cr^, of unit length so (7i = ±1, and the inter- 
action, G, is then equal to G/4. The quenched variable 
ei is equal to 1 if there is a chirality on site i, and is 
equal to if there is no chirality. The lattice sites with 
chirality are to be chosen at random with 50% probabil- 
ity. We impose the additional contraint that the total 
number of chiralities is exactly half the number of lattice 
sites. Defining N to be the number of chiralities we have, 
for every sample. 
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We now discuss the scaling theory and the details of 
our Monte Carlo simulations. We run two independent 
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gcG is defined so that it tends to at high temperatures in 
the thermodynamic limit, and tends to unity as T ^ 
if the ground state is non-degenerate. The chiral glass 
susceptibility should be contrasted with the spin glass 
susceptibility, Xsgj defined by 
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Here a and (3 denote the components of the XY spins in 
Eq. (^) and take values x and y. 

Since the transition in this system is expected to be at 
T = the foJlowing scaling forms are expected for Xcg 
and for (/cJi^: 
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Here the exponent Vcg is the correlation length exponent 
and rjcG is related to the ground state degeneracy — if 
the ground state is unique, one expects rjca =0. At a 
critical point, the data for gcc should be independent of 
size. This is a particularly convenient way of locating the 
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FIG. 1. Results for the Binder moment ratio (?cg, defined 
in Eq. (p^, for different sizes and temperatures. Tfie lines 
are guides to the eye. 



FIG. 2. A log-log plot of results for the chiral glass sus- 
ceptibility XcG, defined in Eq. (|l^), for different sizes and 
temperatures. The lines are guides to the eye. 



transition. The spin glass susceptibility in Eq. (|lj) has 
the finite size scaling form 

Xsa^L^-''^^Xsa{L''''^^T) , (18) 

with exponents Vsg and rjsc which are expected to be 
different from the corresponding chiral glass exponents. 
In this paper, we just focus on the chiral glass critical 
behavior. |— ■ 

We use standard methodalS to ensure equilibration of 
the Monte Carlo simulation. Various quantities are com- 
puted both from overlaps between the two replicas and 
from ^-single replica at different times, see Bhatt and 
YoungtHi for details. Typically, for the largest lattice sizes 
we needed about 100000 Monte Carlo sweeps for equili- 
bration at the lowest temperature. For the averaging 
over disorder we took between 1000 and 10000 samples. 
For most of the data points the statistical errors were 
estimated by averaging over all the samples. However, 
for the largest lattice sizes and lowest temperatures we 
divided up the samples into blocks of a few hundred each 
and calculated the errors from the standard deviations of 
the quantities between different blocks. 

The results for gcc as a function of temperature for 
different sizes are shown in Fig. ^ and the corresponding 
results for Xcg are in Fig. |[ The points at T = are 
obtained by exact enumeration of all the states. 

We also calculated the ground state degeneracy for L 
— 4 and L = 6 and found that there is a small non-zero 
degeneracy, which leads to gcc being slightly below 1 at 
T = 0. However, gcc increases slightly with increasing 
size, so we expect that gca will be unity at T = in the 



thermodynamic limit. One should note here that Ray 
and Mooreo, who worked in the phase representation, 
also found that gcc is not very sensitive to the ground 
state degeneracy, and seems to extrapolate to unity. 

In order to do a finite size scaling analysis for gcc it is 
convenient to incorporate one trivial correction to scal- 
ing. As defined in Eq. (|l2|), g^c tends to 1/iV (rather 
than 0) as T — > oo, which is not completely negligible for 
the sizes studied here. We therefore consider the follow- 
ing quantity 

g'^^^E^S_^^g,,[L^/'^caT)^ (19) 

which does vanish at high temperatures for finite N (and 
tends to unity as T — > if gcc does). 

Fig. H shows our finite size scaling plot for g^^, in- 
cluding data for T < 0.55. Data at higher temperatures 
did not scale well and is presumably not in the scaling 
region. It is clear that good data collapse obtains with 
I'cG = 2.0. We tried scaling with other values of Vcg and 
found that data collapse gets worse both for i^cg > 2.0 
and for Vcg < 2.0. Trying various values for i>cg in this 
way we estimate i>cg — 2.0 ± 0.2 from the data for gcc 
in Fig. |. 

Just at it is convenient to put in the correction to gcG 
in Eq. ( p^ in order that it varies between 1 and as 
the temperature changed, it is also useful to perform a 
similar transformation for Xcg- Since Xcg tends to unity 
as T — > 0, we subtract unity. Furthermore, our zero 
temperature data is consistent with rjcc very close to 
zero so, from now on, we will assume that rjca = and 



3 




FIG. 3. A scaling plot of the data in Fig. |l| according to 
Eq. Data for T < 0.55 has been included. 



that XcG = at T = in the thermodynamic limit. 
Hence we will analyze 
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which varies between 1 and as T increases. 

In Fig. ^ we show a finite size scaling plot for x'cg^ 
including data for T < 0.55. The best fit is for — 1-5: 
somewhat lower than that obtained from gcc- A similar 
differetice is also found in the three-dimensional Ising spin 
glassE^I. 

Combining our exponent estimates from gca and Xcg 
we obtain 



i^cG = 1.8 ±0.3 
T]cG = 0.0 ±0.2 . 



(21) 



These results indicate that the chiralities in the two- 
dimensional XY spin glass order with a correlation length 
exponent i/cg which is different from the spin sJg^s cor- 
relation length exponent, for which earlier workufl found 
UsG —nVi This conclusion is in agreement with earlier 
resultaJ'Q. For the future it would be useful to study 
the i/iree-dimensional XY spin glass in the vortex repre- 
sentation to understand whether there is indeed a finite 
temperatuce chiral-glass phase as found in the work of 
KawamuraQ. 
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FIG. 4. A scaling plot of the data in Fig. |2| according to 
Eq. ( |2c| ) with vcg = 1-5. The data used in the plot is for 
T < 0.55. This plot assumes tjcg = 0, which is reasonable 
since the T = results indicate that rjcc is close to, and 
probably exactly equal to, zero. 
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